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Abstract 

Vacuum expectation value of the surface energy-momentum tensor is evaluated for a 
massive scalar field with general curvature coupling parameter subject to Robin boundary 
CN ' conditions on two parallel branes located on (I?-|- l)-dimensional AdS bulk. The general case 

^«0 ' of different Robin coefficients on separate branes is considered. As an regularization proce- 

^-r , dure the generalized zeta function technique is used, in combination with contour integral 

^— V ' representations. The surface energies on the branes are presented in the form of the sums 

of single brane and second brane-induced parts. For the geometry of a single brane both 
regions, on the left (L-region) and on the right (R-region), of the brane are considered. The 
JL . surface densities for separate L- and R-regions contain pole and finite contributions. For an 

flj ' infinitely thin brane taking these regions together, in odd spatial dimensions the pole parts 

,^ . cancel and the total surface energy is finite. The parts in the surface densities generated by 

the presence of the second brane are finite for all nonzero values of the interbrane separation. 
It is shown that for large distances between the branes the induced surface densities give rise 
/\ ' to an exponentially suppressed cosmological constant on the brane. In the Randall-Sundrum 

braneworld model, for the interbrane distances solving the hierarchy problem between the 
gravitational and electroweak mass scales, the cosmological constant generated on the vis- 
ible brane is of the right order of magnitude with the value suggested by the cosmological 
observations. 

PACS number(s): 03.70. -Fk, ll.lO.Kk, 04.62. -hv 



1 Introduction 

The dynamics of fields on anti-de Sitter (AdS) background possesses remarkable properties and 
much of early interest to this spacetime was motivated by the questions of principal nature 
related to the quantization of fields propagating on curved backgrounds. The presence of the 
both regular and irregular modes and the possibility of interesting causal structure lead to 
many new phenomena. The importance of this theoretical work increased when it has been 
discovered that AdS spacetime generically arises as ground state in extended supergravity and 
in string theories. Further interest in this subject was, and continues to be, generated by the 
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appearance of two models where AdS geometry plays a special role. The first model, the so called 
AdS/CFT correspondence, was motivated by Maldacena 1^ (for a review see [2|) and represents 
a realization of the holographic principle. The AdS/CFT correspondence relates string theories 
or supergravity in the bulk of AdS with a conformal field theory living on its boundary. It 
has many interesting formal and physical facets and provides a powerful tool to investigate 
gauge field theories, in particular QCD. The second model, suggested by Randall and Sundrum 
J2j,[l], is a realization of a braneworld scenario with large extra dimensions. Recently it has 
been realized that the introduction of large extra spatial dimensions may provide a solution to 
the hierarchy problem between the gravitational and electroweak mass scales (HI (for reviews in 
braneworld gravity and cosmology see Refs. [H]-[S]). The main idea to resolve the large hierarchy 
is that the small coupling of four dimensional gravity is generated by the large physical volume of 
extra dimensions. Braneworlds naturally appear in string/M-theory context and provide a novel 
setting for discussing phenomenological and cosmological issues related to extra dimensions. The 
model introduced by Randall and Sundrum is particularly attractive. Their background solution 
consists of two parallel flat branes, one with positive tension and another with negative tension 
embedded in a five dimensional AdS bulk IHJ. The fifth coordinate is compactified on S"^ /Z2, and 
the branes are on the two fixed points. It is assumed that all matter fields are confined on the 
branes and only the gravity propagates freely in the five dimensional bulk. In this model, the 
hierarchy problem is solved if the distance between the branes is about 37 times the AdS radius 
and we live on the negative tension brane. More recently, alternatives to confining particles on 
the brane have been investigated and scenarios with additional bulk fields have been considered 
jni,^ni- Apart from the hierarchy problem it has been also tried to solve the cosmological 
constant problem within the braneworld scenario. The problem of the cosmological constant 
has been considered as the most serious mass hierarchy problem in modern particle physics 
(see, for instance, ^J|) and many attempts addressing this fine-tuning issue can be found in the 
literature. The braneworld theories may give some alternative discussion of the cosmological 
constant (see Refs. [I2l and references therein). The basic new ingredient is that the vacuum 
energy generated by quantum fluctuations of fields living on the brane may not curve the brane 
itself but instead the space transverse to it. 

The investigation of quantum effects in braneworld models is of considerable phenomenologi- 
cal interest, both in particle physics and in cosmology. The braneworld corresponds to a manifold 
with dynamical boundaries and all fields which propagate in the bulk will give Casimir-type con- 
tributions to the vacuum energy (for reviews of the Casimir effect see Refs. ^H1;U1])) and as 
a result to the vacuum forces acting on the branes. In dependence of the type of a field and 
boundary conditions imposed, these forces can either stabilize or destabilize the braneworld. 
In addition, the Casimir energy gives a contribution to both the brane and bulk cosmological 
constants and, hence, has to be taken into account in the self-consistent formulation of the 
braneworld dynamics. Motivated by these, the role of quantum effects in braneworld scenarios 
has received a great deal of attention. For a conformally coupled scalar this effect was initially 
studied in Ref. ^^1 in the context of M-theory, and subsequently in Refs. JEI-IHJ ^^ ^ back- 
ground Randall-Sundrum geometry (for the related heat kernel expansions see Refs. [Sll)- The 
models with dS and AdS branes, and higher dimensional brane models are considered as well 
|28j.|33j-|4Uj. For a conformally cupled bulk scalar the cosmological backreaction of the Casimir 
energy is investigated in Refs. [T5].|28].|33].[IT]-[^. 

In the papers mentioned above, the authors consider mainly the global quantities such as 
the total Casimir energy or conformally invariant fields. The investigation of local physical 
characteristics in the Casimir effect, such as expectation value of the energy-momentum tensor, 
is of considerable interest. In addition to describing the physical structure of the quantum field 
at a given point, the energy-momentum tensor acts as the source in the Einstein equations and 
therefore plays an important role in modelling a self-consistent dynamics involving the gravi- 



tational field. In the case of two parallel branes on AdS background, the vacuum expectation 
value of the bulk energy-momentum tensor for a scalar field with an arbitrary curvature coupling 
is investigated in Refs. |3n],|SJ. In particular, in Ref. [^ the application of the generalized 
Abel-Plana formula |S] to the corresponding mode sums allowed us to extract manifestly the 
parts due to the AdS spacetime without boundaries and to present the boundary induced parts 
in terms of exponentially convergent integrals for the points away the boundaries. The inter- 
action forces between the branes are investigated as well. Depending on the coefficients in the 
boundary conditions, these forces can be either attractive or repulsive. On the background of 
manifolds with boundaries, the physical quantities, in general, will receive both volume and sur- 
face contributions and the surface terms play an important role in various branches of physics. In 
particular, the surface counterterms introduced to renormalize the divergencies in the quasilocal 
definitions of the energy for the gravitational field and in quantum field theory with boundaries 
are of particular interest. For scalar fields with general curvature coupling, in Ref. |1S] it has 
been shown that in the discussion of the relation between the mode sum energy, evaluated as 
the sum of the zero-point energies for each normal mode of frequency, and the volume integral of 
the renormalized energy density for the Robin parallel plates geometry it is necessary to include 
in the energy a surface term concentrated on the boundary (see also the discussion in Refs. 
|14j.|46j). Similar issues for the spherical and cylindrical boundary geometries are discussed in 
Refs. @7j,^8 . An expression for the surface energy-momentum tensor for a scalar field with 
a general curvature coupling parameter in the general case of bulk and boundary geometries is 
derived in Ref. 49 . The purpose of the present paper is to study the vacuum expectation value 
of this tensor for a scalar field obeying Robin boundary conditions on a single and two parallel 
branes of codimension one in (D + l)-dimensional AdS spacetime. In particular, we show that 
in the two-brane set up the surface densities induced on a brane by the presence of the second 
brane are exponentially small for large interbrane distances and give rise to naturally suppressed 
cosmological constant in the brane universe. 

The paper is organized as follows. In the next section we show that the vacuum expectation 
value of the surface energy momentum tensor can be expressed via the expectation values of the 
field square evaluated on the branes. By using the Cauchiy's theorem on residues, an integral 
representations for the related zeta functions for both branes are constructed. They contain 
parts due to a single brane when the second brane is absent, and parts which are induced by 
the presence of the second brane. The latters are finite at the physical point. The analytic 
continuation for single plate contributions is constructed in Section |21 for both regions, on the 
right and on the left of the brane. The surface densities for two-branes geometry are investigated 
in Section ^ Various limiting cases are considered and the cosmological constant induced on 
the brane is estimated. In Section Ewe discuss the balance between the separate parts of the 
vacuum energy and show that they satisfy standard thermodynamic relation. The last section 
contains a summary of the work. 

2 Surface energy-momentum tensor and the generalized zeta 
function 

In this paper we consider a massive scalar field (p{x) on background of a (D + l)-dimensional 
AdS spacetime (AdSo+i) with the line element 

ds'^ = Qikdx'dx^ = e~^''°yri^udx''dx'' - dy^, (2.1) 

and AdS radius given by l/fcz). Here rj^^i, = diag(l, — 1, . . . , — 1) is the metric for the D- 
dimensional Minkowski spacetime, i,k = 0,1,..., -D, and /x, z^ = 0,1, . . . ,D — 1. The corre- 



spending field equation reads 

(^''^ViVfe + m' + CR) ^{x) = 0, (2.2) 

where the symbol Vj is the operator for the covariant derivative associated with the metric gik, 
R = —D{D + l)k'jj is the corresponding Ricci scalar, and C is the curvature coupling parameter. 
For minimally and conformally coupled scalars one has (" = and C = Cc = {D — 1)/4D 
correspondingly. Note that by making a coordinate transformation 

z = e''°y/kD, (2.3) 

metric ()2.ip is written in a manifestly conformally-flat form ds"^ = {kj;)z)^'^7]ikdx^dx^ with x^ = 
z. This is the AdSo+i line element in Poincare coordinates. 

Below we will assume that the field obeys mixed boundary conditions on two parallel infinite 
plane boundaries (branes), located at y = a and y = b, a < b: 

Ay + Bydy^ ^{x) = 0, ?/ = o, 6, (2.4) 

with constant coefficients Ay, By. This type of boundary conditions naturally arises for bulk 
fields in the Randall- Sundrum braneworld due to the Z2 symmetry of the model. The presence of 
boundaries modifies the spectrum for the zero-point fluctuations of the scalar field under consid- 
eration. This leads to the modification of the vacuum expectation values of physical quantities to 
compared with the case without boundaries. For the geometry under consideration, the Wight- 
man function and the corresponding vacuum expectation value of the bulk energy-momentum 
tensor are investigated in Ref. 31, (see also Ref. [2Zj for the case of a conformally invariant 
scalar field). The energy- momentum tensor for a scalar field on manifolds with boundaries in 
addition to the bulk part contains a contribution located on the boundary. For an arbitrary 
smooth boundary dMg with the inward-pointing unit normal vector n , the surface part of the 
energy-momentum tensor is given by the formula [IHI 

T^r'^^ =5{x-dM,)nu (2.5) 

with 

Tik = C^^Kik - (2C - l/2)/iifc(^nV/99, (2.6) 

and the "one-sided" delta-function 5{x;dMs) locates this tensor on dMg. In Eq. 1)2. 6|) . Kif^ = 
h[h^V lUm is the extrinsic curvature tensor of the boundary dMg and hn. = Qik + '^i^fc is the 
corresponding induced metric. Here we are interested in the vacuum expectation values of the 
surface energy-momentum tensor ()2.5() on the branes y = a and y = b. 

Let {(/3q,(x), (/9* (x)} be a complete set of positive and negative frequency solutions to the 
field equation ()2.2|) . obeying the boundary condition 1)2. 4() . Here a is a collective index for all 
quantum numbers. By expanding the field operator over the eigenfunctions ^Pa{x), using the 
standard commutation rules and the definition of the vacuum state, for the vacuum expectation 
value of the surface energy-momentum tensor one finds 



ik 

where 



(0|7;f "^f^lO) = 6{x;dMs){0\nk\0), (2.7) 

{o\Tik\o) = Y,nky^{x),^Ux)}. (2.8) 



Here |0) is the amplitude for the corresponding vacuum state, and the biUnear form Tii.{ip,ip} 
on the right of the second formula is determined by the classical energy- momentum tensor (|2.6() . 
For the geometry of two parallel branes on AdS bulk one has (j = a, b) 

f^ rf (2-9) 

KJJJ = -n^'hng^u, /i, i^ = 0, 1, . . . , L» - 1, 

and Kf)£^ = 0, where n"" and K^'jJ are the inward-pointing unit normal and the extrin- 
sic curvature tensor for the brane at y = j, j = a,b (we consider the region between the 
branes, a < y < b). Note that boundary conditions (|2.4|) can be written in the covariant form 
{Aj + n^^>Bjn^J'^Vi)^p = 0. By using relations (|2.9I) and the boundary conditions, the vacuum 
expectation value of the surface energy-momentum tensor on the brane at y = j is presented in 
the form 

{0\rjfj\0) = -g^.n^^^ [ckD - (2C - l/2)A,/B,] (0|^2|q^^^^^^ (2.10) 

and (0|r|^£)|0) = 0. From the point of view of physics on the brane, Eq. (|2.10jl corresponds to 
the gravitational source of the cosmological constant type, 

(0|rO> |0) = diag (^"^ , -^"^ ,..., -pf^''^) , (2.11) 

with the surface energy density ef^^ (surface energy per unit physical volume on the brane at 



z, 



(surf) 



or brane tension) , stress Pj and the equation of state 



_(surf) (surf) 



ef'^''' = -p'f '''''. (2.12) 

Of course, this is the direct consequence of the Poincare invariance of the branes. 

For an untwisted bulk scalar in the {D + l)-dimensional version of the Randall-Sundrum 
braneworld with brane mass terms Ca and c^, the ratio of the coefficients in the boundary 
condition (|2.4|) is determined by the expression (see, e.g., Refs. |9] . |2l] . [3T] ) 

^ = \ ■ (2.13) 

For the coefficient in formula 1)2. lUj) this gives 

2C - (4C - 1)-% = 8DC(C - Cc) + (4C - 1) n(.')-^. (2.14) 

In particular, the surface energy in the Randall-Sundrum braneworld vanishes for minimally and 
conformally coupled scalar fields with zero brane mass terms. Note that in the supersymmetric 
version of the model ^ one has Ch = —Ca- 

AdS spacetime is divided by the branes into three regions with y < a, a < y < b, and 
y >b. Below in this section, we will consider the region between the branes. The corresponding 
quantities for the other regions are obtained from those as limiting cases and are investigated 
in the next section. Note that in a S^ jZi version of the model, the only bulk is between the 
branes. Due to the L>-dimensional Poincare invariance, the corresponding eigenfunctions can be 
presented in the form 

-<^'^%5Fr' ^"^ <"■'"• (.15) 



UJ 



V /c^ + u^, k = |k|. 



Here the separation constants u are determined by boundary conditions (|2.4I) and will be given 
below. Substituting eigenfunctions (|2.15|) into the field equation H2.2() . one obtains the equation 
for the function fuiu) with the solution 

fu{y) = C„e^'=^J'/2 [j^(^^) ^ b,Y,{uz)] , (2.16) 

in the region a < y <b. Here Jij{x), Yi,{x) are the Bessel and Neumann functions of the order 



{D/2y-D{D + l)C + m2/kl. (2.17) 

The parameter u must be real to ensure stability |Hn],|HJ|. For a given (^ this imposes a lower 
bound for the mass (Breitenlohner-Freedman bound). Note that on AdS bulk the parameter m^ 
can be negative. In the case of a conformally coupled massless scalar, ( = d^c, one has z^ = 1/2 
and the cylinder functions in Eq. H2.16() are expressed via the elementary functions. 

In the region between the branes, the coefficient 5jy is determined from the boundary condition 
on y = a: 

Here and below we use the barred notation 

f(^)(x) = AjF{x) + BjxF'{x), (2.19) 

for a given function F{x), with the coefficients 

Aj = Aj + BjkoD/2, Bj = Bjko- (2.20) 

Note that for a bulk scalar field in the Randall- Sundrum model from Eq. 1)2. 13|) one has 
2Aj/Bj = D{1 — 40 — n^^'Cj/kf). From the boundary condition on the brane y = b we re- 
ceive that the eigenvalues for u have to be solutions to the equation 

gl'^'HuZa^uz,) ^ jl^''\uZa)Y^'\uzb) - Y^''\uZa)Ji'\uzt,) = 0. (2.21) 

This equation gives the spectrum of Kaluza-Klein masses. We denote hy u = Uy^n-, n = 1, 2, . . ., 
the zeros of the function g^ {uZa, uz},) in the right half-plane of the complex variable u, arranged 
in the ascending order, u^,n < iii/,n+i- The set of quantum numbers specifying the eigenfunctions 
is a = (k, n). In our analysis we will assume that the values Aj/Bj are such that there are no 
imaginary zeros. This is the case, for example, when Aa/Ba < and Ah/Bf, > (see also Ref. 
|52j^. The coefficient Cq in Eq. ()2.16|) is determined from the orthonormality condition 

/ dye^^-''^''^yfu^Jy)fu^^^,{y) = 6nn', (2.22) 

Ja 

and is equal to 

^2 _ ^ YJ''\uZa)YJ''\uZf,) _ 

^ "■ -§^9l '[UZa,UZb) 

As it follows from formula 1)2. 1U() . the vacuum expectation values of the surface energy- 
momentum tensor can be obtained from the vacuum expectation values of the field square 
evaluated on the branes. Substituting the eigenfunctions (|2.15p into the corresponding mode 



sum and integrating over the angular part of the vector k, for the expectation value of the field 
square in the region between the branes one finds 

a 

= TTk^-'z^Pp HdhY, ^, ''"''''^'"^ ^="tf (^-'"^i'"--"^) ^ (2.24) 

•^0 n=l Jul n + k"^ -i:9u {uZa, UZb) 

V u=Ujj^n 

where we have introduced the notations 

Pd = ^ n-} ,„ ., , (2.25) 

and 

g^J\u,v) = Mv)%^'\u) - Y,{v)jy\u). (2.26) 

By using the relation gl {u,u) = 2Bj/it, for the corresponding vacuum expectation value on 
the brane at y = j we obtain 

{^W\xm^=^^=2kZ-'zfB,Po rfA:fc^-2V^^ ^ (2.27) 



■^0 n=l Jul + k"^ -iigl {uZa, UZb 



where j,l = a,b, and / ^ j. Quantity (|2.27|) and, hence, the surface energy-momentum tensor 
diverge and need some regularization. Many regularization techniques are available nowadays 
and, depending on the specific physical problem under consideration, one of them may be more 
suitable than the others. In particular, the generalized zeta function method [HSIjE^ is in 
general very powerful to give physical meaning to the divergent quantities. There are several 
examples of the application of this method to the evaluation of the Casimir effect (see, for 
instance, j54j-|57j and references therein). Here we will use the method which is an analog of 
the generalized zeta function approach. 
Instead of (|2.27|) we define the function 

F,{s) = 2k^-hf^ r dkk^'^Y.^uln + ky'' ^.,ng|^(n..n.^n..nz,) (^.28) 

^ -^0 n=l ailS'^ {uZa,UZb)\u=u^,r. 

with ^ an arbitrary mass scale which has been introduced to keep the dimension of the expression. 
Evaluating the integral over k, this expression can be presented in the form 

where B{x,y) is the beta function and we have defined the generalized zeta function as 

oo D+s (0 / \ 

E'^i/,n 9^ \Uv,nZl,Uii^nZj ) 

n=l Q^9i' \UZa,UZb)\u=u„^n 



/. / N V^ ^'^,n 3v [Ui/^nZl,Ui/^nZj) 



We could include in the definition of the zeta function an additional factor n^^^^^"^ to keep 
this function dimensionless. However, this will not affect on the final result for the analytic 
continuation of the function Fj{s). The computation of vacuum expectation value of the surface 
energy-momentum tensor requires the analytic continuation of the function Fj{s) to the value 
s = —1 (here and below |s=_i is understood in the sense of the analytic continuation), 

{Q\^^\Q),=,^=F,{s)\s=-i. (2.31) 



The starting point of our consideration is the representation of the function (|2.3U|1 in terms 
of contour integral: 



Cjis) 



1 



duu 



D+s 



gl'{uzi,uzj) 

(ab) , N 

gl '{uZa,UZb) 



(2.32) 



2TTi J (J 

where C is a closed counterclockwise contour in the complex u plane enclosing all zeros Ui,^n- 
The location of these zeros enables one to deform the contour C into a segment of the imaginary 
axis {—iR, iR) and a semicircle of radius R, R ^ oo, in the right half-plane. We will also assume 
that the origin is avoided by the semicircle Cp in the right half-plane with small radius p. For 
sufficiently large s the integral over the large semicircle in Eq. (|2.H2|) tends to zero in the limit 
R -^ oo, and the expression on the right can be transformed to 



O(^) 



1 

'2Tri 



D+s gi iuzi,uzj) 

' ""^ "" (ab) I N 

Cp gl '{uZa,UZb) 



duu 



1 TT 

+ -sm-{D + l + s) 
vr 2 



duu 



D+s G^j\uZl,UZj) 
Gv '{uZa,UZh) 



(2.33) 



where {j = a, b) 



Gi^\u,v) = l4v)Ki^\u)-K,{v)fJ\u), 
G<^fHu,v) = Ki-Hu)li'Hv) - li''Hu)Ki'\v) 



(2.34) 
(2.35) 



and we have introduced the Bessel modified functions Iu{u) and Kj^{u). Below we will consider 
the limit p ^ 0. In this limit the first integral on the right of formula 1)2. 33() vanishes at the 
physical point s = — 1 , and we will concentrate on the contribution of the second integral. The 
corresponding expression can be presented in the form 

1 TT r°° 

■n>^'Uju{uzj) + BjQjiy{uZa,uzb) 



0(s) ^ 

with the notations 
and 



sm-{D + l + s) 
vr 2 



duu 



D+s 



(2.36) 



Uau{x) 



Ri'^Ux) 



^au{u,v) 



Qbu{u,v) 



Ubu{x) 



K^\v) 



lujx) 

li'hx) 



(2.37) 



.(a) 



(ab). 



Kr{u)Gr>iu,v) 

-li'\v)df\u,v) 



(2.38) 



Now by using the formula 
D-l 



B 



D 



_ Um-{D + l + s) 



vrr(^) 



-f)r(i2^)' 



2 ' 2 J 2- ' T 

the corresponding expressions for the functions Fj{s) can be rewritten in the form 



F,{s) 



'-" -^-'z?B, 



"D 



r(-f)r(i^i^)^ 



s+l 



duu 



D+s 



n^^'Uju{uZj) + BjVljy{uZa, UZb) 



(2.39) 



(2.40) 



-2 ; ^^ J p 

The contribution of the second term in the square brackets is finite at s = — 1 and vanishes 
in the limits Zq ^ and z;, — s- oo. The first term in the square brackets of this expression 
corresponds to the contribution of a single brane at z = Zj when the second brane is absent. 
The regularization is needed for this term only and the corresponding analytic continuation has 
been done in the next section. 



3 Surface energy-momentum tensor for a single brane 

The single brane at z = Za divides the AdS spacetime into two regions corresponding to z < Za 
(L-region) and z > Za (R-region). The properties of the vacuum in these regions are different 
and the corresponding quantities we will differ by the indices R (right) and L (left), respectively. 
Let us consider these regions separately. 

3.1 L-region 

From Eq. ()2.4U|) it follows that for the geometry of a single brane located at z = Za the function 
Fa{s) = Fa (s) is determined by the expression 



l-D 



'^1 r (-}) r (i^^) (^^.)i+" y„ /(°)(.„) 

This integral representation is valid in the strip — (D+1) < Re s < —D and under the assumption 
that the function lif (u) has no real zeros. The latter corresponds to the absence of imaginary 
zeros for the eigenfunction fu{y) in Eq. (|2.15j) with respect to u. This condition is satisfied for 
Aa/Ba > —v- For the analytic continuation to s = — 1, we write the integral on the right of this 
formula as the sum of the integrals over the intervals (0,1) and (1, oo). The first integral is finite 
at s = — 1. To find an analytic continuation of the second integral we employ the asymptotic 
expansions of the Bessel modified function for large values of the argument (see, for instance, 
|58|). For Ba ^ from these expansions one has 



I^{u) 1 ^wi{u 



wE^- M 



where the coefficients wi{i') are combinations of the corresponding coefficients in the expansions 
for the functions Iy{u) and Il{u). The first four coefficients are as follows: 

1 A 
wo{iy) = 1, wi{u) = - - -^, 

2 Ba 

Now we subtract and add to the integrand of the integral over (l,c«) the N leading terms of 
the corresponding asymptotic expansion and exactly integrate the asymptotic part. By this way 
Eq. p.lj) may be written in the form 



(4^)^fcg"^(^z,)-^- f^ l\u.,D^sUn) 



f2^) is) = ^ \ 7 ,r"\ {Ba dun' ,, 



oo 

D+s 



+ I duu 



Iu{u) y-^ Wl{v) 






!i^\u) W'^" 



E wi{u) 
D + s-l 

1=0 



(3.4) 



For N > D — 1 both integrals on the right are finite at s = —1 and Fa (s) has a simple pole at 
s = — 1 corresponding to the summand / = D — 1 in the last sum. As a result, the expression for 
the vacuum expectation value for the field square on the brane contains pole and finite parts: 

{o\^'\om.. = {^'m..,,+{^'t2.^,, (3.5) 



Laurent-expanding the expression on the right of Eq. (|3.4|) near s = —1, one finds 





i^^m... 


= -^k^-' 


for the pole part, and 




i^'tLr- 


Jo 




/•CX) 


' Iu{u) 

''-&\u) 



WD-l{v) 



Pd, 



s + l 



(3.6) 



-1 luju) 



N 

E 

1=0 



+ WD^l{l^) 



\n{txza) + 2^ ( y 



wi{v) 



i*(- 

2^ \2 



AT 

E 

/=0 



wi[y) 



D-l-1 



(3.7) 



for the finite part, with ip{x) being the diagamnia function. By using relations (|2.1U|) and 1)3. 5|) . 
we can write analogous decomposition for the surface energy density in the L-region: 



e 



(surf)(L) _ ^(surf){L) , ^(surf)(L) 



-a.p 



+ e 



a,f 



e 



(surf)(L) 



CkD -{2C- 1/2) Aa/Ba {^') 



o2\(L) 



(3.8) 
(3.9) 



with s = p, f , and for the surface stresses via the equation of state. Note that in the principal 
part prescription adopted in Ref. ^2\ (see also [E^) the Casimir energy corresponds to the 
principal part of the corresponding Laurent expansion near the physical point. In our case this 



is presented by the quantity e^'^'^ 



e 



(surf){R) 
a,f 



for the R-region, see below). 



3.2 R-region 

For a single brane at z = Za the function Fa{s) = Fa (s) in the R-region is determined from 

Eq. (tnni) : 



i^fH^ 



,D-1 



{4TT)^k^-'B^ 



T 



-f)r(i^i^)(/.zji- 



duu 



D+s 






(3.10) 



and differs from the corresponding function in the L-region by the replacement I^ — s- K^, (in 
general, the coefficients in the boundary condition could be different for the left and right 
surfaces of the brane). This integral representation is valid in the strip — (D + 1) < Res < —D 
and under the assumption of absence of real zeros for the function Ku {u) . The latter is the 
case for Aa/Ba < i'. The analytic continuation of the integral on the right to s = —1 can be 
found by the way similar to that for the L-region. For Ba ^ the corresponding asymptotic 
expansion of the subintegrand for large values of the argument has the form 






^ oo 

—Y 



i-iywiiu) 



Ba 



,1+1 



with the same coefficients wi{v) as in Eq. 1)3. 2|) . Now Eq. (|3.1U|) is written in the form 



(3.11) 



i^f)(^ 



{A^T) — k'^-\^^Za)- 



Or( 



D+l+s\ 



B, 



+ 



duu 



D+s 



a I duu 


N 



D+s 



KJu) 



KJf'Uu) 



Ba^TT. + 2^ -ITl 



Ki''\u) 



1=0 



N 



+E 

1=0 



D + s-l 



(3.12) 
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For iV > D — 1 the integrals on the right of this formula are finite at the point s = —1 and 
Fa (s) has a simple pole at s = —1 presented by the summand / = D — 1 in the last sum. 
Hence, the expression for the field square contains pole and finite parts: 

(0|^^|0)iE. = (^^)i^Lp + (v'^)£l,f> (3-13) 

where the separate contributions are obtained by Laurent-expanding expression (|3.12j) near the 
point s = —1: 



(9^^)E.,p = 2(-l)^fc^-^^D-lMf^: 



(3.14) 



and 



N'T / z=Za,i 



2k^-^(3D+i{Ba I duu^' 



/oo 
duu 



D-1 



Ba 





KJu) 



-1 K,{u) 



N 






i-lYwiiu) 



,1+1 



^ j^ i-iYw^ii^) 



1=0 



D-l-1 



(3.15) 



+ {-l)''wD-i{i^) 



ln{fXZa) + -1p i — 



2^ 12 



With the help of these formulae we obtain the similar decomposition for the surface energy 
density in the R-region: 



e 



(surf)(R) _ ^(surf)(R) , -(surf)(R) 



'a,p 



+ e 



a,f 



e 



(surf)(R) 



Ckry-i2C-l/2)AjBa {ip'y^{^,,, 



(3.16) 
(3.17) 



with s = p, f . In the minimal subtraction scheme the pole term appearing in Eq. ()3.16|) is 
omitted. This corresponds to the principal part prescription of Ref. j57j for the total Casimir 
energy. Note that the structure of the pole terms in Eqs. 1)3. 8() and p.l6() allows us to absorb 
them into the corresponding counterterms. There is the freedom to perform finite renormaliza- 
tions since the renormalization counterterms shall be fixed by imposing some renormalization 
conditions (for the detailed discussion of the renormalization procedure for the vacuum energy 
in the Randall-Sundrum braneworld see Refs. PUII^O^ ^. 

The total surface energy density for a single brane at y = a is obtained by summing the 
contributions from the L- and R-regions: 



s. 



(surf)(LR) 



(surf) (L) , (surf) (R) 



(3.18) 



Now comparing the pole parts for the separate regions, we see that in odd spatial dimensions 
these parts cancel out. In particular, taking N = D — 1, for the total surface energy density one 
obtains the formula 



4surf)(LR) =kDf^^^^ l^^j^^ _ (4^ _ 1)^^ 



duu 



D-l 



I,y{u) Ky{u 



/oo 
duu 



D-l 



I,{u) ^ K,{u) 



D-3 

2 



/i-'i..) kI;\u) 



-y 



W2l+l{v) 



Kthu) 



B. 



,21+2 



1=0 



D-3 

2 



-^E 



1=0 



W2l+l{l^) 

D-21-2 



(3.19) 
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where the coefficients wi{v) are defined by relation (|3.2|) . Note that this quantity does not 
depend on the renormahzation scale ^ and the position of the brane. In Fig. ^we have plotted 
the surface energy density (|3.19|) as a function on the ratio Aa/Ba for massless minimally (left 
panel) and conformally (right panel) coupled scalar fields in L) = 3. 





Figure 1: Surface energy density, /c^ Ea [see formula (|3.19l) ]. induced on a single plate as a 
function on the ratio Aa/Ba for massless minimally (left panel) and conformally (right panel) 
coupled scalar fields in D = 3. 



4 Surface densities for two-branes geometry and induced cos- 
mological constant 

In this section we investigate the surface densities generated on a brane by the presence of the 
second brane. These quantities are finite for all nonzero values of the interbrane separation and 
are not affected by finite renormalizations. Below we will concentrate on the energy densities as 
the corresponding stresses on the branes are directly obtained from the equation of state p. 12)1 . 
On the base of relation 1)2. 40() . for the region Za < z < z^ the vacuum expectation value of the 
field square on the brane at y = j (j = a, b) can be presented in the form 



m^\x)\o). 



POO 

Jo 



ji,(UZa,UZb), 



(4.1) 



,(J) 



where {0\(p {x)\0)z=z- is the corresponding quantity for a single brane at y = j when the second 
brane is absent, J = R for j = a and J = L for j = b. The second term on the right of Eq. 
(|4.H1 may also be obtained by using the expression for the corresponding Wightman function 
derived in Ref. ;3lj in the coincidence limit and evaluating its value on the branes. With the 
help of formulae (|2.10j) , (|4.H) , the surface energy density on the brane at z = Zj is presented as 
the sum 

(4.2) 



(surf) (surf)(J) , a (surf) 



where e^^^ is the surface energy density induced on the corresponding surface of a single 



brane at y = j when the second brane is absent (see previous section) , and the term 



Ae 



(surf) 



,(i)( 



2CBj - (4C - 1) Aj n^^\kDZjrBj(3D+i / duu^'-'^j^iuza^uzi,) 





(4.3) 



is the energy density induced by the presence of the second brane. Note that this part is located 
on the surface z = Zq + for the brane at z = Za and on the surface z = Zb — for the brane at 
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z = zi,. On the surfaces 2 = z^ — and 2: = Zf, + the surface densities are the same as for single 
branes. The expression on the right of Eq. H4.3() is finite for all values Za < z\, and is a function 
on the ratio zi,/ Za only. Note that this ratio is related to the interbrane distance by the formula 

V^„ = e'^°(''-'^). (4.4) 

The main contribution into the integral in Eq. (|4.3|1 comes from the modes with u < [zb — Za)~^ 
and the contribution of higher modes is exponentially suppressed. Due to this suppression, the 
same results for the induced densities will be obrained in the model where instead of externally 
imposed boundary condition the fluctuating field is coupled to a smooth background potential 
that implements the boundary condition in a certain limit |59j . 

Using the Wronskian for the Bessel modified functions, it can be seen that 



[S|(x2z2 + i,2 



A]] %.l 



•^^a} ^'^bj 



nS-'' Zj— — In 



-&\ 



XZa)Kl \xZh) 



rW/ 



Iv'{xZb)Kl '[XZa 



.(a). 



(4.5) 



for j = a,b. This allows us to write the expression (|4.,3j) for the surface energy density in another 
equivalent form: 



. (surf) _i,DD+ln a 



duu 



^„i 2QB, + (1 - 40^,- 



_5_ 

dzi 



In 



ir {uZa)Kr {UZ,) 



/O -y 



42 



#. 



UZb)Kif' (uZa) 



(4.6) 



This form of the surface energy density will be used below in the discussion of the energy balance. 
Let us consider the limiting cases of the part (|4.3() . For large values of AdS radius to 
compared with the interbrane distance, ku{b — a) <C 1, the main contribution to the integral 
on the right of Eq. ()4.3() comes from the large values of uza ^ [koib — a)]~^. Assuming that 
Ba/{b — a) and m{b — a) are fixed we see that the order of the Bessel modified functions is large. 
Replacing these functions by their uniform asymptotic expansions for large values of the order 
(see (HH])) one obtains 



Aep^«2n(^)(l-4C)i,-4/?D+i 



du- 



u^iu^ 



m 



D_ 

I 2 ' 



HA] 



Byr^ 



{Aa-Bau){Ai,+B,,u) ^2u(b-a) _ i 

{Aa+BaU){Ai-Bi,u) 



(4.7) 



The expression on the right is the corresponding surface energy for two parallel Robin plates 
in the Minkowski spacetime. It can be easily checked that for a massless scalar field the sum 



Aea""" + Ae^'^'^'^ evaluated from Eq. (|4.7|) coincides with the result obtained in Ref. ^]. In 
the limit of small interbrane separation, (b — a) -^ 0, for a fixed k£), by using the asymptotic 
formulae for the Bessel modified functions for large values of the argument, from the general 
formula in D > 2 dimensions we find 



^^(surf) 



2C - (4C - l)Aj/Bj 



n 



U) 



(TjkDCRiD - i)r 



'D+l^ 



2D-H^{D 



l){b 



\D-1 



where Cr{^) is the Riemann zeta function, aj = 1 for |i?a/^a|, |^6/^b| ^ ^d(& — 
1 for \Bj/Aj\ ^ kj:)(b — a) and Bi/Ai = 0, with I = b for j = a and / 



i2-D 



(4.8) 

, and 
a for 



j = b. We see that for small interbrane distances the sign of the induced surface energy density 

is determined by the factor 2C, — (4^ — l)Aj/Bj and this sign is different for two cases of aj. 

Now we turn to the asymptotic behavior of the surface energy density 1)4. 3|) for large distances 
between the branes, kr){b — a) » 1. Introducing a new integration variable x = uzb, by making 
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use of formulae for the Bessel modified functions for small values of the argument, and assuming 
\^a\ 7^ |-Sa|i^, to the leading order we receive 



.D+2uAk^Ba 



2CBa - (4C - 1) Aa 



A.^^^- 1^ : 1 ::, ' ^..--^^-^%^. (4.9) 



2z^) T,^V{%)V'^{y){Aa-Bavf }^ Il'\x) 



n+2u-i ^i^ {^) 



and 



2u 



Aei-'^f) = fcg ( ^ ) f^#, (4.10) 



where we have introduced the notation 

B, \2CB, - (4C - 1) it 

J 1/ 



00 ^D+2u-l 



rf^^TMT— (4-11) 



22u+D-i^iY {§) uF^u) Jo Ti^^\x) 

For \Aa\ = \Ba\i^ it is necessary to take into account the next terms in the corresponding 
expansions of the Bessel modified functions for small values of the argument. The integral in 
Eq. ()4.9() is negative for small values of the ratio Ai,/Bi, and is positive for large values of 
this ratio. As it follows from Eq. 1)4. 10(1 . the sign of the quantity AeJ,''^'^ for large inter btane 
separations is determined by the combination {Blu"^ — A^) [2^ + (1 — ^QAi^/Bi)] of the coefficients 
in the boundary conditions. 

For large values of the mass with m ^ ki) and m{b — a) ^ 1, as we see from Eq. ()2.17() . 
one has u ~ m/kD ^ 1- Introducing a new integration variable u = vy and using the uniform 
asymptotic expansions for the Bessel modified functions, to the leading order from (|4.3|) one 
receives 

Ae(-rf) ~ ^(,) 2Ci?, - (4C - l)i, {Ai+n^')mBi){mkD)''l^mB^ ^-2Mb-a) u 12) 

where I = a for j = b and / = b for j = a, and we keep mBi an arbitrary. As we could 
expect in this case the induced surface densities are exponentially suppressed. In Fig. |21we have 
plotted the surface energy densities determined by Eq. 1)4. 3|1 as functions on Za/z^ and Ah/Bi, 
for Aa/Ba = -2. 

To discuss the physics from the point of view of a D-dimensional observer residing on the 
brane y = j, it is convenient to introduce rescaled coordinates x'-^ in accordance with x'-'^ = 
Q-^DJx^. With these coordinates the warp factor in the metric is equal to 1 at the brane y = j 
and they are physical coordinates for an observer on this brane. Now after the dimensional 
reduction of the action, by the way similar to that in the Randall-Sundrum braneworld (see, 
e.g., 0)) it can be seen that Z)-dimensional Newton's constant Gdj measured by an observer on 
the brane at y = j is related to the fundamental (D + l)-dimensional Newton's constant Gd+i 
by the formula (for the gravitational equation on the brane see, e.g., |H].[5n|) 

Gdj = 4^^^?^l^%±^e(^-2)^-(^-^). (4.13) 

Note that in the orbifolded version of the model an additional factor 2 appears in the denom- 
inator of the expression on the right. For large interbrane distances one has Gua ~ k^Go^i, 
Gob '^ k£)GD+ie>^~^'^°^^~°'\ and the gravitational interactions on the brane y = b are ex- 
ponentially suppressed. This feature is used in the Randall-Sundrum model to address the 
hierarchy problem. Now we will show that this mechanism also allows to obtain a naturally 
small cosmological constant generated by vacuum quantum fluctuations. 
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Figure 2: The surface energy densities {z^/ Zj)^ ^e^^' /k^, j = a,b, as functions on Za/z;, and 
Ab/Bb for a massless minimally coupled scalar field with Aa/Ba = —2 in D = 3 induced on the 
branes at z = Zq (left panel) and z = zi, (right panel). 



The surface energy density (|4.3j) corresponds to the gravitational source of the cosmological 
constant type induced on the brane at z = Zj by the presence of the second brane. Within two 
brane geometry discussed in this section, from the point of view of an observer living on the 



brane at z = Zj 
by the relation 



the effective cosmological constant induced by the second brane is determined 



Anj = SvtGdj Ae 



(surf) 



SvrAe 



(surf) 



M 



D-2 
Dj 



(4.14) 



where Mdj is the Z)-dimensional effective Planck mass scale for an observer on the brane. By 
using relation (|4.13)) . for the numerical estimates it is convenient to present the ratio of the 



induced cosmological constant (|4.14l) to the corresponding Planck scale quantity in the brane 
universe in the form 



A 



Dj 



Ae 



(surf) 



SttGdjM^^ 






Zb 



D 



k 



D 



Mda 



o(o-i) 

D-2 



D 



{Zh/Za) 



D-2 



1 



D 
D-2 



(4.15) 



rl-D 



where M/)+i is the fundamental {D + l)-dimensional Planck mass, Gd+i = Mj^,^ . For large 
interbrane distances, by taking into account formulae (|4.9|) . (|4.10j) . one obtains the following 
estimate 



A 



Dj 



D+2u 



zJ \Md+i) 



D(D-l) 
2 



(4.16) 



SnGDjMg^ 

showing that this ratio is of the same order of magnitude for both branes. However, due to 
the large hierarchy between the Newton's constants on the branes, for the ratio of the induced 
cosmological constants one has hj^^/^Da ~ {zb/ Za)"^- 

In the Randall-Sundrum {D + l)-dimensional braneworld the brane at z = z^ corresponds 
to the visible brane. For large interbrane distances by taking into account Eq. (|4.10|) . for the 
ratio of the induced cosmological constant (|4.14|1 to the corresponding Planck scale quantity in 
the corresponding brane universe one obtains 



Ad6 


\Zb, 


D+2v 


( ko ^ 


D{D-1) 
0-2 BaV + Aa 


^T^GDbM^," 


\Md+i^ 


BaV - Aa 



fib) 



(4.17) 
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where the function fj) is defined by Eq. (|4.11|) . As this function determines the value of the 
induced cosmological constant for large interbrane distances, in Figure 01 we have plotted its 
dependence on the ratio of the coefficients in the boundary condition on the visible brane and 
on the mass of the scalar field for minimally and conformally coupled cases. 



0.2 



-0.1 





0.06 



0.04 



0.02 



Figure 3: The quantity /^ defined by Eq. (|4.11|) as a function on the ratio Ah/Bb and m/kD 
in the model with D = 4 for minimally (left panel) and conformally (right panel) coupled scalar 
fields. 

The higher dimensional Planck mass M^^i and AdS inverse radius k^ are two fundamental 
energy scales in the {D + l)-dimensional AdS bulk space which in the Randall-Sundrum model 
are usually assumed to be of the same order, ku ~ M^+i (see, e.g., ;6_). In this case from Eq. 
(|4.17|) one obtains the cosmological constant which is exponentially suppressed to compared with 
the corresponding Planck scale quantity on the visible brane. In the original Randall-Sundrum 
braneworld with D = A, Md+i ~ TeV, Mub = Mpi ~ 10^^ TeV, to account for the observed 
hierarchy between the gravitational and electroweak scales, in accordance with Eq. ()4.13|1 we 
need Zb/za ~ 10^^ (by formula (|4.4j) this corresponds to the interbrane distance about 37 times 
larger than the AdS radius k]^ ). In the case of a scalar field with the mass |m^| < kj^ (recall 
that on AdS bulk the parameter m^ can also be negative), for the induced cosmological constant 
this leads to the estimate which is of the right order of magnitude with the value implied by the 
cosmological observations. 

5 Energy balance 

In this section we will consider the balance between the separate parts of the vacuum energy. 
In the region between the branes the total vacuum energy per unit coordinate volume on the 
branes is the sum of zero-point energies of elementary oscillators: 



2 



d^-'k 



(2vr) 



D- 



l / j\k + Uiy^n, 



1/2 



(5.1) 



n=l 



where u^^n are roots of equation (|2.2H) . The total vacuum energy in the bulk region Za < z < Zb 
is obtained by the integration of Q-component of the volume energy-momentum tensor over this 
region: 

i?(™i)= /ci^+ixx^(0|ri™')°|0). (5.2) 
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Using the formula for the volume energy- momentum tensor from Ref. |31| . we can see that this 
energy differs from (|5.1() (see Ref. |49| for the discussion in general case of bulk and boundary 
geometries). This difference is due to the presence of the surface energy located on the branes. 
The surface energy per unit coordinate volume on the brane, E^^'^^^' , is related to the surface 
densities from Eq. (|2.11|) by the formula 

(surf) 
J=a,b ■' 

Now it can be easily checked that the total vacuum energy is the sum of the volume part and 
the surface part: 

Expression (|5.1|) may be evaluated using the procedure which is basically the same as that 
used in Refs. ^J (see |551 I56j for reviews) for studying the Casimir effects in the geometry of 
spherical boundaries. Within the framework of the Randall-Sundrum braneworld, this has been 
done in Refs. |18l I2L)1 I24j by the dimensional regularization method and in Ref. :22; by the zeta 
function technique. Refs. |18| I2fl| |22' consider the case of a minimally coupled scalar field in 
Z? = 4, and the case of arbitrary C, and D with zero mass terms c^ and c^ is calculated in Ref. 
|24j . Here we briefly outline the zeta function approach in the general case. Instead of (|5.1|) we 
introduce the following zeta function 

where, as above, the parameter ^ with dimension of mass is introduced by dimensional reasons. 
Evaluating the integral over k one receives 

,,s+l p ( s-D+l \ oo 



C(-) = ^\n^ V^ ' E<n'-'- (5.6) 



(47r)^ r^l 



n=l 



We need to perform the analytic continuation of the sum on the right of this formula to the 
neighborhood of s = — 1. Transforming it into a contour integral and by deforming the contour 
appropriately, the total energy in the region Za < z < Zf, is presented in the form 

E = las)\s=-i = E%^ + E%^ + Ai?, (5.7) 

where E^^za ^^ ^^^ vacuum energy in the region z > Za for a single brane at z = Za, and E),J^^ 
is the vacuum energy in the region z < Zb for a single brane at z = Zb- The interference term 
AE in this formula is finite for all nonzero values of the interbrane separation (see below, Eq. 
(|5.1(i|) ) and the analytic continuation is needed for single brane parts only. Under the assumption 
Bb ^ 0, for the region z < z^ the latter is given by the expression 



'^"' 2r (f ) r {^^) Jo du 



^^''h 



UZb) 



(5.8) 

s=-l 



where we have introduced the notation 



^^J''\uzb) = J^^I('\uzb), (5.9) 

UZb Bb 
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and, as before, |s=_i is understood in the sense of the analytic continuation. The expression for 
the energy EJ^^ is obtained from Eq. H5.9|) by the replacement 

9 pUZa 



Sp)(nz,)^S(^''^)(nz,) 



TTUZa Ba 



■Ki''\uZa) 



After the analytic continuation we can see that E)^^ contains pole and finite parts. 
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1=1 



D-l 



+ Dw^Siu) 



1=1 



ln(/iz;,) + ^v(y + l 



(5.13) 



1 

-i' 

2^ 



(i). 



for the finite part, with N > D — 1. In these formulae, wi (u) are the coefficients in the 



asymptotic expansion of the function In Sj> ' (u) for large values of the argument: 



In 



oo 

1=1 



W 



ij)l 



a, b. 



(5.14) 



These coefficients can be related to the coefficients in the similar expansions for the Bessel 



modified functions (see, for instance. 



The corresponding formulae for the energy E_°^^ are 

obtained from expressions (|5.12j) and (|5.13j) by replacements (|5.10j) and w^ (u) -^ (— l)'ii;[" (z^). 
The renormalization of the divergences in the corresponding formulae for the vacuum energies 
can be performed by using the brane tensions (for a detailed discussion of the renormalization 
procedure within the framework of the Randall-Sundrum model see, e.g., Refs. HO] , jH] , [ISl ) • 

As in the case of the surface energies, now we see that in the calculation of the total vacuum 
energy for a single brane at z = Za in odd spatial dimensions, including the contributions from 
L- and R-regions, the pole parts of the energies cancel out (assuming that the coefficients in the 
boundary conditions (|2.4|) on the right and left surfaces are the same) and we obtain a finite 
result. In particular, taking N = D — 1, for this energy one receives 
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(5.15) 
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The expression on the right of this formula can be easily evaluated numerically. Note that the 
energy per unit physical volume on the brane is determined as {kDZa)^E^"'' and does not depend 
on the brane position. 

Unlike to the single brane parts, the interference term AE in Eq. 1)5. 7|) is finite and is 
determined by the formula: 
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(5.16) 



Note that this part of the vacuum energy is not affected by finite renormalizations. Now let us 
check that this quantity obeys the standard energy balance equation. We expect that in the 
presence of the surface energy this equation will be in the form 



dE = -pdV+ Y,^t''^dS^'\ 



j=a,b 

where V is the [D + l)-volume in the bulk and S^^' is the D-volume on the brane y 
coordinate volume on the brane, 

rb 
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dye-^^^\ S^^^ 



-Dkoj 



J = a, b, 



(5.17) 



j per unit 



(5.18) 



In Eq. 1)5. 17() . p is the perpendicular vacuum stress on the brane and is determined by the vacuum 



expectation value of the |)-component of the bulk energy-momentum tensor: p 
Combining equations 1)5. 17() and (|5.18|) . one obtains 
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with pi^' being the perpendicular vacuum stress on the brane at y = j. In Ref. [49_' it has been 
shown that these stresses can be presented as sums of the self-action and interaction parts: 
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j = a,b. (5.20) 

j when the second brane 



The first term on the right is the pressure for a single brane at y 

(i) 
is absent, and p\^^) is induced by the presence of the second brane. The latter determines the 

interaction forces between the branes and is defined by the formula |49j 
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(5.21) 



Assuming that relation 1)5. 19(1 is satisfied for the single brane parts, for the interference parts 
one obtains 
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(5.22) 



Now by taking into account expressions ()4.6I) , (|5.16l) , (|5.21|) for the separate terms in this formula, 
we see that this relation indeed takes place. Hence, we have checked that the vacuum energies 
and effective pressures on the branes obey the standard energy balance equation. Note that here 
the role of the surface energy is crucial and the vacuum forces acting on the brane [determined 
by P/fnt)]' ™ general, can not be evaluated by a simple differentiation of the total vacuum energy. 
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6 Conclusion 

The natural appearance of AdS in a variety of situations has stimulated considerable interest 
in the behavior of quantum fields propagating in this background. In the present paper we 
have investigated the expectation value of the surface energy-momentum tensor induced by the 
vacuum fluctuations of a bulk scalar field with an arbitrary curvature coupling parameter sat- 
isfying Robin boundary conditions on two parallel branes in AdS spacetime. The Wightman 
function and the vacuum expectation value of the bulk energy-momentum tensor for this ge- 
ometry are investigated in Ref. |HJ (see also Ref. 1^0] for the case of the Randall-Sundrum 
braneworld). By making use the expression for the surface energy- momentum tensor from Ref. 
|49j and the boundary conditions on the branes, the expectation values of the surface energy 
density and stresses are expressed via the expectation values of the field square on the branes. 
As an regularization procedure for the latters we use the generalized zeta function technique, in 
combination with contour integral representations. Using the Cauchy's theorem on residues, we 
have constructed an integral representations for the zeta functions on both branes, which are 
well suited for the analytic continuation. These functions are presented as sums of two terms. 
The first ones correspond to the zeta functions for single branes when the second brane is absent. 
The second terms are induced by the presence of the second brane and are finite at the physical 
point. For the analytic continuation of the single brane zeta functions we subtract and add 
to the integrands leading terms of the corresponding asymptotic expansions, and present them 
as sums of two parts. The first one is convergent at the physical point and can be evaluated 
numerically. In the second, asymptotic part the pole contributions are given explicitly. As a 
consequence, the single brane surface Casimir energies for separate L- and R-regions contain 
pole and finite contributions. The remained pole term is a characteristic feature for the zeta 
function regularization method and has been found in the calculations of the total Casimir en- 
ergy for many cases of boundary geometries. As in the case of the total vacuum energy, the 
renormalization of these terms can be performed by using the brane tensions. For an infinitely 
thin brane taking L- and R-regions together, in odd spatial dimensions the pole parts cancel 
and the surface Casimir energy is finite. In this case the total surface energy per unit physical 
volume on the brane (surface tension) does not depend on the brane position and can be directly 
evaluated by making use formula 1)3. 19(1 . The results of the corresponding numerical evaluation 
for minimally and conformally coupled scalar fields are presented in Figure ^ As seen from 
the corresponding graphs, in dependence of ratio of coefficients in the boundary condition the 
surface energy for a single brane can be either negative or positive. The cancellation of the pole 
terms coming from oppositely oriented faces of infinitely thin smooth boundaries takes place 
in very many situations encountered in the literature. It is a consequence of the fact that the 
second fundamental forms are equal and opposite on the two faces of each boundary. In even 
dimensions there is no such a cancellation. 

The surface densities induced on the brane by the presence of the second brane are inves- 
tigated in Section [l] The corresponding contributions to the zeta functions are finite at the 
physical point giving finite quantities for all nonzero values of the interbrane distance and are 
not affected by finite renormalizations. The corresponding energy densities are determined by 
formula (|4.3j) and are located on the surfaces corresponding to the region between the branes. 
They can also be presented in another form given by Eq. (|4.6() . In the limit of large AdS 
radius. A;/) — > 0, from this formula the result for two parallel plates on the Minkowski bulk is 
recovered. For large distances between the branes, the surface densities induced by the second 
brane are exponentially suppressed by the factor exp[/c£)(2i/ + D){a — b)] for the brane at y = a 
and by the factor ex.p[2iyk£){a — b)] for the brane at y = b. The exponential suppression also 
takes place in the large mass limit. From the viewpoint of an observer living on the brane at 
y = j, the D-dimensional Newton's constant is related to the higher-dimensional fundamental 
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Newton's constant by formula H4.13|l and for large interbrane separations is exponentially small 
on the brane y = b. The corresponding effective cosmological constant generated by the second 
brane is determined by Eq. (|4.14|) and is suppressed to compared with the corresponding Planck 
scale quantity in the brane universe by the factor exp[k£i{2i' + D){a — 6)], assuming that the 
AdS inverse radius and the fundamental Planck mass are of the same order. In the original 
Randall-Sundrum model with D = 4, for a scalar field with the mass \rin?\ < k'^, and interbrane 
distances solving the hierarchy problem, the value of the cosmological constant on the visible 
brane by order of magnitude is in agreement with the value suggested by current cosmological 
observations without an additional fine tuning of the parameters. 

The vacuum energy localized on the boundary plays an important role in the consideration 
of the energy balance. In Section [3 we consider the total vacuum energy in the region between 
the branes, evaluated as a sum of zero-point energies for elementary oscillators. It is argued that 
this energy differs from the energy, obtained by the integration of the bulk energy density over 
the region between the branes. We show that this difference is due to the presence of the surface 
energy located on the branes. Further we briefly outline the procedure for the regularization of 
the total vacuum energy by the zeta function technique. This energy is presented as a sum of 
single branes and interference parts. The single brane terms contain pole and finite contributions 
and explicit formulae are given for both these parts. In the calculations of the total vacuum 
energy for a single brane in odd spatial dimensions the pole parts of the energies, coming from 
L- and R-regions, cancel out and a finite result emerges. This energy per unit physical volume 
on the brane is independent of the brane position. For the geometry of two parallel branes, 
the interference part in the total vacuum energy is given by formula (|5.16|) and is finite for all 
nonzero values of the interbrane separation. Further, we have shown that the induced vacuum 
densities and vacuum effective pressures on the branes satisfy the energy balance equation 1)5. 17() 
with the inclusion of the surface terms, which can also be written in the form (|5.22jl . 

Acknowledgments 

The author acknowledges the hospitality of the Abdus Salam International Centre for Theoretical 
Physics, Trieste, Italy. This work was supported in part by the Armenian Ministry of Education 
and Science, Grant No. 0887. 

References 

[1] J.M. Maldacena, Adv. Theor. Math. Phys. 2, 231, (1998); Int. J. Theor. Phys. 38, 1113 
(1999). 

[2] O. Aharony, S.S. Gubser, J. Maldacena, H. Ooguri, and Y. Oz, Phys. Rep. 323, 183 (2000). 

[3] L. Randah and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999). 

[4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999). 

[5] N. Arkani-Hamed, S. Dimopoulos, and G. Dvah, Phys. Lett. B 429, 263 (1998); Phys. Rev. 
D 59, 086004 (1999); I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. Dvah, Phys. 
Lett. B 436, 257 (1998). 

[6] V.A. Rubakov, Phys. Usp. 44, 871 (2001). 

[7] P. Brax and C. Van de Bruck, Class. Quantum Grav. 20, R201 (2003). 

[8] R. Maartens, |gr-qc/Q312Q59| 

21 



[9] T. Gherghetta and A. Pomarol, Nucl. Phys. B586, 141 (2000). 

[10] W.D. Goldberger and M.B. Wise, Phys. Rev. D 60, 107505 (1999); A. Pomarol, Phys. Lett. 
B 486, 153 (2000); H. Davoudiasl, J.L. Hewett, and T.G. Rizzo, Phys. Lett. B 473, 43 
(2000); Y. Grossman and M. Neubert, Phys. Lett. B 474, 361 (2000); R. Kitano, Phys. Lett. 
B 481, 39 (2000); C.V. Chang and J.N. Ng, Phys. Lett. B 488, 390 (2000); R. Altendorfer, 
J. Bagger, and D. Nemeschansky, Phys. Rev. D 63, 125025 (2001); S.R. Huber and Q. Shah, 
Phys. Rev. D 63, 045010 (2001). 

[11] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989); V. Sahni and A. Starobinsky, Int. J. Mod. 
Phys. D 9, 373 (2000); S.M. Carroll, Living Rev. Rel. 4, 1 (2001); N. Straumann, 
|astro-ph/0203330; P.J.E. Peebles and B. Ratra, Rev. Mod. Phys. 75, 559 (2003); T. Pad- 
manabhan, Phys. Rept. 380, 235 (2003). 

[12] N. Arkani-Hamed, S. Dimopoulos, N. Kaloper, and R. Sundrum, Phys. Lett. B 480, 
193 (2000); S. Kachru, M. Shulz, and E. Silverstein, Phys. Rev. D 62 045021 (2000); 
C.P. Burgess, R.C. Myers, and F. Quevedo, Phys. Lett. B 495, 384, (2000); S. F5rste, 
Z. Lalak, S. Lavignac, and H.P. Nilles, Phys. Lett. B 481, 360 (2000); S. F5rste, 
Z. Lalak, S. Lavignac, and H.P. Nilles, J. High Energy Phys. 09, 034 (2000); C. Shmid- 
huber, Nucl. Phys. B585,385 (2000); J.-W. Chen, M.A. Luty, and E. Ponton, J. High 
Energy Phys. 09, 012 (2000); P. Binetruy, J.M. Cline, and C. Grojean, Phys. Lett. 
B 489, 403 (2000); S.P. de Alwis, Nucl. Phys. B597 263 (2001); C. Csaki, J. Elrich, 
and C. Grojean, Nucl. Phys. B604, 312 (2001); N. Tetradis, Phys. Lett. B 509 307 
(2001); S.-H.H. Ty and I. Wasserman, Phys. Rev. Lett. 86, 1682 (2001); E. Flanagan. 
N. Jones, H. Stoica, S.-H.H. Ty and I. Wasserman, Phys. Rev. D 64, 045007 (2001); 
J.E. Kim, B. Kyae, and H.M. Lee, Phys. Rev. Lett. 86, 4223 (2001); D.K. Park, H. Kim, 
S. Tamaryan, Phys. Lett. B 535, 5 (2002); C.L. Gardner, Phys. Lett. B 524, 21 (2002); 
G. Dvah, G. Gabadadse, and M. Shifman, Phys. Rev. D 67, 044020 (2003); S. Nojiri 
and S.D. Odintsov, JCAP 06, 004 (2003); S. Khlebnikov, Phys. Lett. B 551, 211 (2003); 
U. Giinter, P. Moniz, and A . Zhuk, Phys. Rev. D 6 8, 044010 (2003); S.M. Carroll and 
M.M. Guica, hep-th/0302067; M. Setare, hep-th/0308109 ; Li-Xin Li, lhep-th/0310075! 
G. Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, and S. Zerbini, 'hep^th/0312269 
Y. Aghababaie, C.P. Burgess, S.L. Parameswaran, and F. Quevedo, Nucl. Phys. B680, 



389, (2004); C.P. Burgess, Eep-th/0402200 ; J. Vinet and J.M. Cline, hep-th/0406141 



[13] V.M. Mostepanenko and N.N. Trunov, The Casimir Effect and Its Applications (Oxford 
University Press, Oxford, 1997); G. Plunien, B. Muller, and W. Greiner, Phys. Rep. 134, 
87 (1986); The Casimir Effect. 50 Years Later, edited by M. Bordag (World Scientific, 
Singapore, 1999); M. Bordag, U. Mohideen, and V.M. Mostepanenko, Phys. Rep. 353, 
1 (2001); M. Bordag, ed.. Proceedings of the Fifth Workshop on Quantum Field Theory 
under the Influence of External Conditions, Int. J. Mod. Phys. A 17 (2002), No. 6&7; 
K.A. Milton, The Casimir Effect: Physical Manifestation of Zero-Point Energy (World 
Scientific, Singapore, 2002). 



[14] K.A. Milton, |hep-t h/0406"024l 



[15] M. Fabinger and P. Horava, Nucl. Phys. B580, 243 (2000). 

[16] S. Nojiri, S. Odintsov, and S. Zerbini, Phys. Rev. D 62, 064006 (2000). 

[17] S. Nojiri and S. Odintsov, Phys. Lett. B 484, 119 (2000). 

[18] D.J. Toms, Phys. Lett. B 484, 149 (2000). 

22 



[19] S. Nojiri, O. Obregon, and S. Odintsov, Phys. Rev. D 62, 104003 (2000). 

[20] W. Goldberger and I. Rothstein, Phys. Lett. B 491, 339 (2000). 

[21] S. Nojiri and S. Odintsov, J. High Energy Phys. 07, 049 (2000). 

[22] J. Garriga, O. Pujolas, and T. Tanaka, Nucl. Phys. B605, 192 (2001). 

[23] A. Flachi and D.J. Toms, Nucl. Phys. B599, 305 (2001). 

[24] A. Flachi and D.J. Toms, Nucl. Phys. B610, 144 (2001). 

[25] I. Brevik, K.A. Milton, S. Nojiri, and S.D. Odintsov, Nucl. Phys. B599, 305 (2001). 

[26] A. Flachi, I.G. Moss, and D.J. Toms, Phys. Lett. B 518, 153 (2001); Phys. Rev. D 64, 
105029 (2001). 

[27] A.A. Saharian and M.R. Setare, Phys. Lett. B 552, 119 (2003). 

[28] E. Elizalde, S. Nojiri, S.D. Odintsov, and S. Ogushi, Phys. Rev. D 67, 063515 (2003). 

[29] J. Garriga and A. Pomarol, Phys. Lett. B 560, 91 (2003). 

[30] A. Knapman and D.J. Toms, Phys. Rev. D 69, 044023 (2004). 

[31] A.A. Saharian, prep-th/0312092 



[32] M. Bordag and D.V. Vassilevich, J. Phys. A 32, 8247 (1999); LG. Moss, Phys. Lett. B 491, 
203 (2000); P.B. Gilkey, K. Kirsten, and D.V. Vassilevich, Nucl. Phys. B601, 125 (2001). 

[33] S. Nojiri, S. Odintsov, and S. Zerbini, Class. Quantum Grav. 17, 4855 (2000). 

[34] N. Kaloper, J. March-Russel, G.D. Starkman, and M. Trodden, Phys. Rev. Lett. 85, 928 
(2000); G.D. Starkman, D. Stojkovic, and M. Trodden, Phys. Rev. Lett. 87, 231303 (2001); 
G.D. Starkman, D. Stojkovic, and M. Trodden, Phys. Rev. D 63, 103511 (2001). 

[35] W. Naylor and M. Sasaki, Phys. Lett. B 542, 289 (2002). 

[36] I.G. Moss, W. Naylor, W. Santiago- German, and M. Sasaki, Phys. Rev. D 67 125010 (2003). 

[37] A. Flachi and O. Pujolas, Phys. Rev. D 68, 025023 (2003). 

[38] A. Flachi, J. Garriga, O. Pujolas, and T. Tanaka, J. High Energy Phys. 08, 053 (2003). 

[39] A.A. Saharian and M.R. Setare, Phys. Lett. B 584, 306 (2004). 



[40] J.P. Norman, hep-th/0403298 



[41] S. Mukohyama, Phys. Rev. D 63, 044008 (2001). 

[42] R. Hofmann, P. Kanti, and M. Pospelov, Phys. Rev. D 63, 124020 (2001). 

[43] A.H. Yeranyan and A.A. Saharian, Astrophysics 46, 386 (2003). 

[44] A.A. Saharian, "The generalized Abel-Plana formula. Applications to Bessel functions and 



Casimir effect," Report No. IC/2000/14, |hep-th /0002239 
[45] A. Romeo and A.A. Saharian, J. Phys. A 35, 1297 (2002). 
[46] S.A. Fulling, J. Phys. A 36, 6857 (2003). 



23 



[47] A. A. Saharian, Phys. Rev. D 63, 125007 (2001). 

[48] A. Romeo and A. A. Saharian, Phys. Rev. D 63, 105019 (2001). 

[49] A. A. Saharian, Phys. Rev. D 69, 085005 (2004). 

[50] P. Breitenlohner and D.Z. Freedman, Phys. Lett. B 115, (1982); Ann. Phys. (N.Y.) 144, 
249 (1982). 

[51] L. Mezincescu and P.K. Townsend, Ann.Phys. (N.Y.) 160, 406 (1985). 

[52] A. A. Saharian and M.R. Setare, Class. Quantum Grav. 20, 2765 (2003). 

[53] J.S. Dowker and R. Critchley, Phys. Rev. D 13, 3224 (1976); S.W. Hawking, Commun. 
Math. Phys. 55, 133 (1977); G.W. Gibbons, Phys. Lett. B 60, 385 (1977). 

[54] E. Ehzalde, S.D. Odintsov, A. Romeo, A. A. Bytsenko, and S. Zerbini, Zeta Regularization 
Techniques with Applications (World Scientific, Singapore, 1994). 

[55] K. Kirsten, Spectral Functions in Mathematics and Physics (CRC Press, Boca Raton, FL, 
2001). 

[56] A. A. Bytsenko, G. Cognola, E. Elizalde, V. Moretti, and S. Zerbini, Analytic Aspects of 
Quantum Fields (World Scientific, Singapore, 2003). 

[57] S.K. Blau, M. Visser, and A. Wipf, Nucl. Phys. B310, 163 (1988). 

[58] Handbook of Mathematical Functions, edited by M. Abramowitz and LA. Stegun (Dover, 
New York, 1972). 

[59] N. Graham, R.L. Jaffe, V. Khemani, M. Quandt, M. Scandurra, and H. Weigel, Phys. Lett. 
B 572, 196 (2003); N. Graham, R.L. Jaffe, V. Khemani, M. Quandt, M. Scandurra, and 
H. Weigel, Nucl. Phys. B645, 49, (2002). 

[60] T. Shiromizu, K. Maeda, and M. Sasaki, Phys. Rev. D 62, 024012 (2000); C. Barcelo and 
M. Visser, J. High Energy Phys. 10, 019 (2000). 

[61] E. Elizalde, S. Leseduarte, and A. Romeo, J. Phys. A 26, 2409 (1993); S. Leseduarte and 
A. Romeo, J. Phys. A 27, 2483 (1994); A. Romeo, Phys. Rev. D 52, 7308 (1995); M. Bordag, 
E. Ehzalde, and K. Kirsten, J. Math. Phys. 37, 895 (1996); M. Bordag, K. Kirsten, and 
J.S. Dowker, Commun. Math. Phys. 182, 371 (1996); M. Bordag, E. Elizalde, K. Kirsten, 
and S. Leseduarte, Phys. Rev. D 56, 4896 (1997). 



24 



